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I. INTRODUCTION 

The use of vertex functions derived from the covariant Bethe-Salpeter (BS) equation [jl],0] 
or the covariant spectator (or Gross) equation [0,H to calculate matrix elements involving 
bound states requires that the normalization of these vertex functions be determined. Al- 
though this may seem like a trivial problem, the complicated forms of these equations, along 
with the generally nonlocal nature of their interaction kernels makes it more than trivial. 
In fact we know of only one derivation of the normalization condition for the three-body BS 
bound state vertex functions [§], and this derivation does not seem to show that the normal- 
ization condition follows simply from application of the wave equation to the description of 
scattering in the region of the bound state, and does not obtain the normalization condition 
for the Faddeev subvertex functions. The three body spectator equations have only recently 
been derived in their final form [^] and this paper includes the first derivation of the nor- 
malization condition of the three body spectator vertex functions. Normalization conditions 
for the two-body BS and Gross equations have been previously obtained by many people 
(see, e.g., 0-0), but this paper includes a first demonstration that charge is automatically 
conserved by the spectator equations. 

Following this introductory discussion the paper begins, in Sec. II, with a derivation of 
the normalization condition for both the two-body BS and Gross vertex functions. These 
results, which are not new, are presented in order to introduce the techniques which will 
be used later in the derivation of the three-body conditions. At the conclusion of this 
section we prove, for the first time, that the normalization of the two-body Gross vertex 
function is identical to the requirement that the charge of the bound state, as defined by the 
charge operator introduced by Gross and Riska |]TU[, is equal to the sum of the charges of 
its constituents. This shows that charge is automatically conserved by the spectator theory. 

The third section contains the new derivations of the normalization conditions for the 
three-body Bethe-Salpeter and Gross vertex functions [Q, which closely follow the approach 
used in the two-body case. A discussion of the Feynman rules used in this paper and the 
conventions used in symmetrizing the scattering matrices is presented in the Appendix. 

II. NORMALIZATION OF RELATIVISTIC TWO-BODY VERTEX FUNCTIONS 

The derivations of the normalization condition for two and three-body vertex functions 
are very similar. In this section we start with the simpler two-body case, and then generalize 
the derivation to the three body case in Sec. III. 

A. Normalization of the BS two-body vertex function 

The two-body Bethe-Salpeter equation for the scattering matrix Ai is represented by 
the Feynman diagrams in Fig. 0. Using the Feynman rules described in the Appendix, this 
corresponds to 

M = V -VGbsM (2.1) 

= V-MGbsV, (2.2) 
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FIG. 1. Diagrammatic representation of the two-body BS equation for the scattering matrix. 



where V is the two-body interaction kernel and Gbs is the free two-body propagator. In 
terms of the one body propagator Gi defined in Eq. (|2.3CI| ) below, the BS propagator is 

GBs = -tGiG2. (2.3) 

From Eq. ( p.2| ) we have 

V = M+MGbsV, (2.4) 

and substituting this equation into Eq. (pTl| ) gives the following nonlinear equation for TW 

M = V- MGbsM - MGbsVGbsM . (2.5) 

If the two-body system has a bound state at P^ = M^, where P is the total four- 
momentum, the scattering amplitude Ai will have a pole at P^ = M^. Therefore, the 
two-body scattering amplitude can be written 



M 



|r) (r| 

M2 - P2 



+ 7^, 



(2.6) 



where |F) is the Bethe-Salpeter vertex function and (F| is the representation of |F) in the 
dual space, including Dirac conjugation. The first term is the contribution from the bound 
state propagator, with possible spin degrees of freedom associated with the propagation of 
the bound state suppressed, and 7^ is a finite at P^ = M^. Substituting this form ( |2.6| ) into 
the linear form (p.l|) of the equation for A^, and equating residues at P^ 
wave equation for the vertex function |F) 

|F) = -VGbs\T). 



M gives the 



(2.7) 

Substituting (|2.6| ) into the nonlinear form of the scattering equation ( p.5|) , and keeping 
all terms which are singular near the bound state pole at P^ = M^, gives the following 
relation 



|r) (r| 



lim 



|r) 



{T\GBsil + VGBs)\T) 



(r| 



A'P - P2 

-UGbs {l + V Gbs) \T) {T\ - |F) (F|(l + Gbs V) Gbs 7^} • 



(2.8) 



Near the bound state pole at P^ = M ^ the BS equation ( p.TD ensures that the last two terms 
involving the remainder TZ vanish. Therefore, as P^ — > M^, Eq. (^.8|) reduces to 



Since both the numerator and the denominator vanish in this hmit, we expand the numerator 
around P^ = M^, giving 

1 = -(r| -^ |r) - (r| ^ (GssVGss) |r) . (2.10) 

(The BS equation ensures that any terms proportional to d\r)/dP'^ will also vanish.) Dis- 
tributing the derivative in the second term over the product Gbs^Gbs and using the bound- 
state BS equation ( p.7|) reduces the normalization condition to 

1 = (r| G'j,s |r) - (r| GbsV'Gbs |r) . (2.11) 

where 

, ^dGB^^P^dGBS_ d^^P^dV_ 

^Bs Qp2 2p2 5p/. ' ap2 2P2 apM ' ^ ■ ^ 

and all derivatives are taken at the bound state pole Pq = E{P) = (P^ + M^)^/^. 

For identical particles, the symmetrized scattering matrix, which we will denote by M, 
must satisfy 

V12M = MV12 = (M (2.13) 

where P12 is the permutation operator that exchanges all labels for particles 1 and 2, and 
(^ = 1 for bosons and ( = —1 for fermions. This can be achieved by symmetrizing the 
unsymmetrized scattering matrix A4 according to 

M = A2M = A2MA2 , (2.14) 

where 

^2 = ^(1 + Cn2) (2.15) 

is the symmetrization operator and the last relation follows from A2M. = M.A2 and A^ = 
A2- Using the Bethe-Salpeter equation for the scattering matrix yields 

M = A2M =A2{V- VGbsM) 

= A2V - A2VGBSA2M 

= A2V-A2VGbsM (2.16) 

where the identities A2V = A2VA2 and A2GBS = Gbs-A~2 have been used in rewriting the 
expression. Defining the symmetrized kernel illustrated diagrammatically in Fig. ^ 

V = A2V=]^{l + CVi2)V, (2.17) 
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FIG. 2. Diagrammatic representation of the symmetrized kernel. 

gives a Bethe-Salpeter equation for the symmetrized scattering matrix with the same form 
as the Bethe-Salpeter equation for nonidentical particles 



M = V~ VGbsM 



(2.18) 



Since the structure of this equation is the same as in the case of identical particles, the 
derivation of the bound state equation and the normalization condition is the same as before 
with replacement oiV ^V 



where 



s{T\G'Bs\n-s{T\GBsV'GBs\n 



|r)s = A|r) 



(2.19) 



(2.20) 



is the symmetrized vertex function. 

Note that the symmetrized scattering matrix defined above is normalized differently from 
the symmetrized scattering matrices presented in most field theory texts, as described in 
the Appendix. 



B. Normalization of the Gross t^vo-body vertex functions 

The Gross or spectator equation for distinguishable particles can be obtained from the 
corresponding Bethe-Salpeter equation by means of a simple prescription. For unequal mass 
particles with relatively long-range interactions, the heavier of the two particles is placed on 
its positive-energy mass shell in performing the energy loop integral over intermediate state 
four- momenta. Here, the heavier particle is assumed to particle 1. The scattering matrix 
then becomes 



M = V -VQ1G2M 



(2.21) 



where the projection operator Qi = Qj places particle i on its positive-energy mass shell. 



Equation (|2.21|) does not represent a closed set of solvable equations. In order to obtain 
these it is necessary to also place particle 1 on shell in both the initial and final states leading 
to 



where Va = QiVQi and Mu = QiMQi- This is illustrated in Fig. H. 



(2.22) 
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FIG. 3. Diagrammatic representation of the two-body Gross equation for the scattering matrix. 
The X on the hne for particle 1 indicates that it is on shell. 



The situation for identical particles is somewhat more complicated than for the Bethe- 
Salpeter equation. Since ^12^2^12 = Gi and 'Pi2Qi'Pi2 = Q2, the scattering matrix can no 
longer be symmetrized by simply symmetrizing the final state. The solution to this problem 
is to use a symmetric version of the spectator equation 



M = V-lv (Q1G2 + Q2G1) M 



(2.23) 



where the factor of | is necessary to prevent double counting of the elastic unitary cut of the 
scattering matrix. The factor of | emerges automatically if we think of obtaining Eq. ( p^.23 ) 
from the BS equation by doing the integration over the relative energy variable by averaging 
the results obtained from closing the contour in the upper half and lower half planes (and 
keeping the positive energy nucleon poles only). Eq. ( |2.23|) can now be symmetrized in the 
same way as the Bethe-Salpeter equation, giving 



M = V-- F (Q1G2 + 22^1) M 

= V- ^F (Q1G2 + P12Q1G2P12) M 

= V-^VQiG2il + CVi2)M 
= V -VQ1G2M . 



(2.24) 



Clearly, Q1G2 could be eliminated in favor of 22^*1 in exactly the same manner. Therefore, 
the M-matrix is independent of the choice of on-shell particle, and the resulting integral 
equation for the scattering matrix for ( |2.24| ) can be taken to be 



Mil = V^ii - ViiG2M_ 



11 



(2.25) 



V 



1 

2 



^f^ 



^<r 




FIG. 4. Diagrammatic representation of the symmetrized kernel for Gross equation. As in 
Fig. 3, the x indicates that the particle is on shell. 
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FIG. 5. Diagrammatic representation of the normalization condition for the two-body Gross 
vertex function. The dashed Hne represents the derivative d/dP"^, and the x means that the 
particle is on shell. 

where the symmetrized potential is shown in Fig. ^. 

As for the Bethe-Salpeter equation, the existence of a bound state implies a pole in the 
scattering matrix, 



M- 



11 



iri) {T^ 



M2-P2 



+ 7^. 



(2.26) 



where |Fi) is the bound state vertex function with particle 1 on shell. Repeating the devel- 
opment used in the preceding section gives the Gross equation for |Fi) 



\T,) = -VnG2\Ti). 



(2.27) 



We will choose the momentum of the on-shell particle and the total momentum P as in- 
dependent variables so that d/dP^ refers to differentiation with respect to P^ holding the 
momentum of the on-shell particle constant. Hence the on-shell projection operators Q will 
be independent of P„, and the normalization condition becomes: 



l = (Fi|G"2|ri)-(Fi|G'2V;iG'2|ri) 



(2.28) 



This equation is represented diagrammatically in Fig. ^. Note that this result could also 
be obtained by applying the spectator prescription directly to the unsymmetrized version 
of the Bethe-Salpeter normalization condition ( |2.19| ). Clearly, a similar expression can be 
derived for r2. 



C. Conservation of electromagnetic charge for the Gross equation 



Since we have obtained the normalization condition ( p.28| ) directly from the two-body 
equation, with no reference to the charge of the bound state, we now are in a position to 
prove that the charge of the bound state is conserved; ie. that the charge of the bound state 
is the sum of the charges of its two constituents, regardless of its structure. 

The construction of the matrix elements for the two-body current for the Gross equation 



have been discussed by Gross and Riska [jTO[ and applied to elastic electron scattering from 
deuteron in 0. As argued in 0, the bound state matrix elements of the electromagnetic 
current in the simple case when the OBE kernel is constructed from the the exchange of 
neutral bosons is represented by the Feynman diagrams shown in Fig. |^. Notice that since 
in our spectator formulation only the first particle is being placed on-shell, the first diagram 
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FIG. 6. Diagrams which contribute to the electromagnetic charge of the two body bound state. 



in Fig. 1^ with 1 <-h> 2 does not appear. Here we show explicitly that the last two diagrams 
in Fig. ^ give a contribution in the limit g — >■ equal to that of the first diagram as well as 
the derivative term of the potential as it appears in the normalization condition ( |2.28|) . 



Our argument will require that we replace the operator form of the Gross equation (|2.27| ) 
by its explicit representation in momentum space. We will use the notation 



N 



1 for bosons 
2m for fermions 



(2.29) 



and 



Qi 



A. 



1 for bosons 
for fermions 



Gi{p) 



mf — p^ — ie 



1 



2mj 



TTii— i) — ie 



for bosons 



for fermions . 



(2.30) 



where Aj is what remains of the projection operator Qj after the integral over the relative 
energy has been carried out and the delta function which places particle 1 on shell has 
been removed, and p denotes an on-shell four- momentum with p^ = rn? . The equation then 
becomes 



r(pi;P) 



d^kxN — 



i^-Kf^Ek, 



V{p,, h; P)A,{h)G2{P - h)T{h; P) . 



(2.31) 



Note that the vertex functions are written as a function of the four-momentum of particle 
1, pi, and the total four-momentum, P. Multiplying both sides of the equation by the 
projection operator Ai(pi) and using A^(A;i) = Ai(A;i) gives 



ri(pi;P) 



d^kiN ^ 



{27r)^2Ek, 



Vn{pi,h;P)G2{P-h)T,{h;P) 



(2.32) 



where Ai(pi)V^(pi, ki, P)AiCki) = V^ii(pi, k; P) and Tiipu P) = Ai(pi)r(pi; P). 

Using the Feynman rules, and the notation used in Eq. ( |2.31| ) and ( |2.32| ), the diagrams 
in Fig. E yield 



J^{P\ P) = j .^0^T,{p,; n G,{P' - P,)f,{P' -PuP- Pi)G,{P - pi) ri(pi; P) 
+ Jj0^T{p, + q;P')G,{p, + q)j^{p, + q,p,)G2{P-pi)T,{p,;P) 

+ J J0^^^(P'^-^ P') MPi^Pi - l)Gi{p[ - q)G,{P' - p[) T{p[ ~q;P), (2.33) 

where Ji{p'iP) is the renormahzed e.m. current of i-th constituent particle. When taken 
between the on-shell states for p' = p, this current defines a physical charge Cj. The one-body 
currents satisfy the one-body Ward-Takahashi identity 

%Jt{p'^,P^) = -e. {G7\p') - Gr\p)) ^ -e. q, ^G-\p) , (2.34) 

where G is a full one-particle propagator and the last form of the identity holds near the 
point 5'^ -^ 0. Comparing the terms linear in q^ we get from (|2.34| ) 



3^{p,p) = -ei^GT\p). (2.35) 

Any purely transverse parts Jxiip'-iP) of the off-shell one-body current, which satisfy 
qfiixiip' iP) = and cannot be determined from ( 2.34| ), are assumed to vanish in the limit 



g ^ 0. Therefore, they do not affect the charge e,. Next, note the identity 

G^{p)3^{p,p)G,{p) = -e,G,{p) (J-Gr\p)\ G,{p) = e. ^G,{p) . (2.36) 

Thus, under the assumptions discussed above the insertion of the g = photon into the line 
of the constituent particle is unambigously determined by its physical charge and derivative 
of its propagator. In this paper we always approximate G by (|2.30| ) with a physical mass 
mi] our arguments are however valid also for the generalized Gi{p) = h'^(jp')/{mi— p) used 
in a dynamical model of Ref . [§] . 

Now, the charge of the bound system is found by taking the limit of ( ^.33| ) g — *> 0, or 
P' -^ P, since contributions from any higher multipoles, such as magnetic dipole or charge 
quadrupole must vanish for q = 0. In the limit considered the propagators for particle 1 in 
the last two terms, Gi{pi + q) and Giip'^ — q), become singular. However, it is easy to see 
that these singularities cancel to give a finite result. 



Before we proceed with the evaluation of the g — > limit of Eq. ( 2.33 ), it is instructive to 
describe this cancellation of singularities in very general terms. To this end we note that the 
last two terms of Eq. ( |2.33| ) can be identically represented by the four dimensional integral 

J^iP', P)| ,a. . = -z / j^Tip^ + g; P')G^{p, + q)j^ip, + q,p,) 

terms JC (27rj^ 

xG,{pi)G2{P-pi)T{pi;P), (2.37) 

where the integration over pio is done along the contour G in the complex pio plane, shown 
in Fig. 1^. This contour surrounds the positive energy poles of the propagators Gi (pi -|- g) (at 
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FIG. 7. Location in the pio complex plane of the 6 poles from the three nucleon propagators 
in Eq. (2.34). The last two terms in Eq. (2.33) emerge if the integral of pio over the contour C 
enclosing the poles 3a and la is evaluated using the residue theorem. 

3a) and Gi{pi) (at la), and evaluating the integral by the residue theorem gives only the 
last two terms in Eq. ( |2.33| ), where for notational simplicity in the third term the integral 
over d^pi is replaced by the integral over d^p[ (where p'^^ = Pi + q)- It has been shown |Tl| 
that for elastic reactions the cuts of T{pi; P) never overlap with poles 3a and la of Fig. |^, 
this makes the representation ( 2.37|) possible and unambiguous for all values of momentum 
transfer q. It is now clear that as g ^ the poles at 3a and la coalesce into a double pole, 
giving a finite result. Calculation of the residue of the double pole requires a calculation of 
the derivative of the rest of the integrand, which also explains the appearance of derivatives 
in the final result below. 

We now return to the reduction of Eq. ( p.33|) . Substituting ( p.36|) into the first term of 
Eq. ( |2.33D gives 

J^{P,P)\n.. =e, f ,/f^^jI T,{p,; P') ^G2{P - Pi) ri(pi; P) , (2.38) 



rirst 
term 



pi 



9P„ 



where d/d{P — p\)^i was replaced by d/dP^, which is possible only because pi is an on-shell 
four-momentum, so that all four of its components are independent of P. We will return to 
this expression after we have reduced the last two terms. 

The last two (singular) terms in Eq. ( p.33| ) can be reduced by expressing the vertex 
function with both particles off-shell in terms of the vertex function with particle 1 on-shell 
by using the equation 



r(pi;P) 



^3t. TU 

' -Vip,,h;P)G2iP-h)T,ih;P), 



{27:f2Ek, 



(2.39) 



where the kernel V{pi,ki]P) now connects an incoming channel with particle 1 on-shell 
to an outgoing channel with both particles off-shell and is obtained from V{pi,ki\P) by 
replacing the four momentum pi = {Ep, p) by the off-shell four-momentum (Pq ~ Ep^p, p). 
[Eq. (|2.39|) is obtained from the original Eq. ( p.21|) by going to the bound state pole in the 
usual way.] Using this equation to iterate each of the off-shell vertex functions once yields 



J^^(P',P)|la.t 



2 
terms 



d^p[ N 



d^PiN 



{2Txf2Ep,J (27r)3 2Ep/ 



-Ti {p\- P')G,{P' - p[)j^M, P';Pi. P) 

xG2{P-pi)Ti{pi;P) (2.40) 
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where 



ji:M,p';pi.p) 



£kiN 



-Ai(i^i: 



Vip[,h + q;P')G,Cki + q) 



xj^Ch + q, k)Aiik)G2{P - kWCkuPi; P) + 
V{p[, h; P')G2{P' - k)Aiik)j^Cku k - q) 



xGi{k-q)V{k-q,Pi]P) 



Ai(pi) 



(2.41) 



Using the WT identity (|2.34|) for finite q, and recaUing that Gi^{k)^i{ki) = 0, gives for 
the divergence of the effective current 



q^j^M,P';PuP) = -eiAi{Pi)J 



(fkiN 



{2n)^2Ek, 

x\v{p[,k + q;P')G2{P-k)Mk)V{k,Pi;P) 
-F(p;, k; P')Ai{k)G2{P' - kWik - q,Pi; P) 



Ai(pi) . (2.42) 



Note that this step removes the propagators Gi{ki + q) and 6*1(^1 — q) which are singular 
in the g — »• hmit. The purely transverse contributions to the effective current ( |2.41| ) arise 
only from the transverse parts of the one-body current and so, by assumption, vanish as 
g — i> 0. Therefore, we can use the divergence of the effective current to extract its finite part 
for P' = P. As before we expand ( [2.42 ) about the point q = and equate terms linear in 
q^ to obtain 



j^M,P;pi.P) = -eil 



(PkiN 

{2ny2Ek, 



Ai(p;) 



X 



d 



l\dk 



-Vip[, k;P)] G2{P - k)Aiik)vCki,Pi; P) 



iM 



-V{p\, k; P) {4pG2{P - k)\ Ai{k)V{k,Pi; P) 

^ Vi,{k,Pi;P) 



+Vip[,k;P)G2iP-k)Aiik 



Bk 



iM 



Ai(pi). (2.43) 



Substituting this into Eq. (p.40|) , and using K\ = Ai and the spectator equation ( p.32|) for 
the bound state gives 



d^p^N - 



^''(^•^)l:™' =»'/(5fer.fe;P)U^G.(P-pOJr,fe;P) 



+ei 



Vi[p^-P)G2[P -Pi) ——V{p^,pi;P) 



+ J^^(^i'^i; P)) G2{P - Pi)Ti{pi; P) . 



dp 



iM 



(2.44) 
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Note that the first term is identical to Eq. ( p. 381) if ei ^ 62- To reduce the last two terms 
use the fact that the symmetrized one boson exchange (OBE) potential has the form 



V{p[,PuP) = lil + C'Pl2)V{p[,p^;P) = ^[v,{p^-p[) + CVSl+p'l-P) 



(2.45) 



where Vd is the direct term and Ve the exchange term as defined in Ref. P]. Therefore only 
the derivatives of the exchange term contribute to the charge operator, and 



d 



+ 



d 



A 



5pi^ dp[^J 



V{p[,PuP) = lc^ ^ 



+ 



d 



'l^^ 



VSi+p'i-P) 



'if^. 



-C J^K(pi +P[ -P) = -2^Vip[,p,; P) , 



(2.46) 



where we again used the fact that pi and p[ are both independent of P. Combining all of 
these results together gives our final result for the charge operator of the bound state 



:r^(p,p) = (61 + 62) 



TiiPi; P) — ^75 — ri(pi; p) 



-2ei 



(27r)3 2Ep, 
d^p[d^Pi N^ 
{2Txf AEp.Ep^- 



dP, 

T,{p[; P)G2ip[; P)^^^-^%P^G2{pi; P)T^ip^■, P) 



dp. 



CB^Pr 



(2.47) 



where we used A? = Ai and Ai(p[)[dV{p[,pi; P)/dP^]Ai{pi) = dVu{p[,pi; P)/dP^. In 
the last step we also assumed, for definiteness, that the bound state has spin zero and total 
charge cb- Multiplying the equation by P^ and using the identity 



P., 



dP„ 



2P2 



gp2 



(2.48) 



gives the relation 

es = (ei + 62 

- 261 



d^PiN 
(27r)3 2Ep, 
d^p[d^pi N'^ 



{2'KfAE^>E, 



T,ip,;P)G',ip,;P)T,{p,;P) 

-T,{p[; P)G2{p[; P)v[M,Pu P)G2{pi; P)Upu P) 



pi 



(2.49) 



We now distinguish two cases. If the two particles are not identical, the OBE kernel does 
not include an exchange term and the V^ term vanishes [recall Eq. ( p.46| )]. In this case the 
normalization condition gives 65 = 61 + 62. If the particles are identical, the V^ term will 
not vanish, but the charges will be equal, so 61 + 62 = 26i = cb and we obtain the same 
result. For either identical or distinguishable particles, the normalization condition ensures 
the conservation of charge. 

This proof can be generalized to kernels which include multiple meson exchange. In this 
case the additional energy dependence contained in the kernel will lead to new contributions 
to the V term in the normalization condition, but these new terms will be reproduced 
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exactly by the q ^ limit of additional Feyninan diagrams containing interaction currents 
also required by the multiple meson exchange kernel, and the charge will still be conserved. 
The charge conservation and its relation to the normalization condition can be conve- 
niently discussed in terms of N-particle Ward-Takahashi identities [0]. We will present 



these identities together with the general e.m. currents for two- and three-body systems for 



the spectator formalism in a forthcoming paper |13 



III. NORMALIZATION OF RELATIVISTIC THREE-BODY VERTEX 

FUNCTIONS 

In this final section we find the normalization condition for relativistic three-body vertex 
functions. The derivation for the three-body Bethe-Salpeter equation is formally identical 
to that given above for the two-body case. We briefly repeat it just to set a stage for 
discussion of more complicated spectator case. In particular, we discuss in some detail 
the symmetrization necessary for the identical particles and at the end we express our 
normalization condition in terms of the symmetrized Faddeev subamplitudes. 

A. Normalization of three-body Bethe-Salpeter vertex functions 

1. Distinguishable Particles 
The three-body Bethe-Salpeter equation for distinguishable particles can be written 

r = v-vG%sr (3.1) 

= V-TG%sV, (3.2) 

where G%g = —G1G2G3 is the three body BS propagator, and 

V = z ^ V'Gi^ = E ^' (3-3) 

i i 

is the sum of all products of two body potentials V^ multiplied by the inverse of the prop- 
agator —iGi of the third, non-interacting spectator, and we are assuming that there are 
no relativistic three body forces (ie. diagrams which are three-body irreducible). In these 
equations and those which follow we have adopted the odd-man-out or spectator notation. 
If i, j and k represent some permutation of 1, 2 and 3, the two-body potential V"* represents 
the interaction of particles j and k. 

These equations look deceptively simple, but the singularities in V make them difficult to 
solve, and they are usually reduced to Faddeev form by introducing the subamplitudes T** 
in which particle i' is the initial spectator and particle i the final spectator. Since these ob- 
jects represent components of a summed Feynman perturbation expansion of the scattering 
problem, the terms "initial" and "final" refer to the topological form of the corresponding 
Feynman diagrams, not to temporal priority. Dropping the convention of summation over 
repeated indices, so that all sums will he written explicitly, the initial equations for these 
subamplitudes follow immediately from ( |3.1| ) and 
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T''' = i Sn,V'G-^ - V'G'bs E ^^'' (3-4) 

e 

= 1 5u^V'Gi' - J2 r^'GisV , (3.5) 

e 

where the singularities have been canceled using iG^^G^^g = —iGjGk = G^^g. Note that 

VG'^s = V'G^Bs ■ (3-6) 

We will adopt the convention that the index i can take on any value, and that i', j', and 
k' are all different indices describing initial state particles (just as i, j, and k describe final 
state particles). Note that 

T = Y,r'' (3.7) 



and that Eq. ( |3Tl| ) follows by summing ( |3^ ) over i and i', and that Eq. ( p.2| ) follows by 



summing ( p.5|) over i and i' . Convergence of these equations can be further improved by 
moving the T" from the rhs to the Ihs and introducing the two-body scattering amplitude 

M' ={l + V'G'bs)'^ V' = V^^ (l + V'G'bsY^ ■ (^-^^ 
This gives the famous Faddeev equations for the subamplitudes 

T'"' = i 5ii,M'G-^ - M'G'ss E ^^'' (3-9) 

= i du'M'GT^ - E T'^G%sM'' . (3. 10) 



The Feynman diagrams representing Eq. ( p.9|) are shown in Fig. |g. In these figures the solid 
dot denotes the spectator. 

Since the logic behind the resummation of the multiple scattering series in terms of 
Faddeev amplitudes is the same as in the nonrelativistic three-body problem, it is not sur- 



prising that the equations have a form similar to the nonrelativistic Faddeev equations [jT4 
It is important to remember, however, that in this case the summation is done in terms of 
Feynman perturbative expansion rather than the time-ordered perturbation theory of the 
nonrelativistic case. As a result global three-body propagators are replaced by products of 
one-body Feynman propagators and an inverse propagator must be included in the driving 
term to ensure that no propagator appears for external legs of the three-body T matrix. 
Similarly, due to similarities in the underlying logic, many of the expressions derived below 
have nonrelativistic analogues of similar form. The connection between the relativistic for- 
mulation of the three-body problem and the usual nonrelativistic approach will be explored 
elsewhere in detail. 

A bound state of the three-body system is associated with the occurrence of a pole at 
P^ = M ^ in the three-body scattering amplitude T (and in each of the subamplitudes T** ) , 
where P is the total momentum of the three-body system and M is the mass of the bound 
state. Therefore, T and T** can be written as the sum of a pole term and a regular part i?** 

14 



T = 









FIG. 8. Diagrammatic representation of the Faddeev equations for the amphtudes T^*. Note 
that the spectator is identified by the sohd dot. 



T 



|P) (P 



M2-P2 

|r) (r| 






(3.11) 
(3.12) 



where |r) = J2i |r*). Inserting the ansatz ( p.ll|) into the scattering equations and taking 
the limit P^ -^ M^ gives the equation for the bound state vertex subamplitudes 






(3.13) 



Similarly, using either Eq. ( p.8| ) or starting with the equations for T, we obtain 

\T) = -Y: V'GU^) = -VG^sl^) . (3.14) 

e 

The derivation of the three-body normalization condition is most easily obtained from 
the ansatz (|3.12| ) and the original equations (|3.1|) and ([31^ ). Repeating the steps which lead 
to Eq. ( ^.11[ ) we obtain for the three-body vertex function the normalization condition of 
exactly the same form 






-(r|G° 



d 



VG° 



BS 



|r). 



(3.15) 



^p2 I / \ ^ Bb Qp2 

The two forms of the normalization condition ( 3.15| ) are identical because of the bound state 
Eq. (§3D. 



2. Identical particles 

The equations and normalization condition for identical particles can be obtained im- 
mediately by symmetrizing the results given above. However, to lay the foundation for the 
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discussion of the spectator equations it is convenient to describe this symmetrization in 
detaiL Most of the operator algebra below is exactly the same as in non-relativistic case. 
The normalized three-body antisymmetrization operator is 

^3 = - (1 + CP12 + C^13 + C^23 + Vi+ Pg) 



1 + CVij) (1 + CV^k + CVjk) 



1 

3! 

^^{l + CV,k + CVju){l + CV,j) , (3.16) 



where 



"^4 — 'Plz'Pl2 — "^23^13 — 'Pl2'P2'i 

V5 = V23VU = V12VU = ^13^23 • (3. 17) 

Note that ^3 = ^3^2 = -4.2^3 where we may choose to have the two body symmetrization 
operator A2 act on any of the three two-body subspaces. We will exploit this property in 
the following discussion. 

Using the relations VjkV^Vjk = V and VjiV^Vij = V^ (which also hold for Gi) we can 
readily show that the symmetrization operator ^3 commutes with the total potential 



The symmetrized version of Eq. (|3.1| ) is then obtained by multiplying the equation by ^3, 
which gives 

T „ 

= AsV - AsVG^sAsT 

(3.18) 

where 



(3.19) 

A similar argument works for the transposed equation ( |3.2D , and hence the symmetrized 
amplitude T satisfies the equations satisfied by T with V replaced by V. Note that we have 
written the kernel (|3.19|) in terms of the symmetrized two-body kernels (|2.17| ), and that V 
can also be written in terms of a single V 

V = ^3 E ^' = ^3 E 'Pi^^''P^i = ^SV ' E C'P^e . (3.20) 

e e e 

As before, i can be any index, and consistency requires the definition Va = C^ so that 
(Vii = 1. The last step in {^^ follows from Vj.V'Vij = V^ and AsVu = (A3. The 



AsT 


= AsV 


- A^Gls'T 




= A3V 


-A^G'^j.sA^T 




= V - VGlsT , 


V = 


--A3V = 


A.Y^V 


r = 


--%Tg: 


-1 



symmetric version of Eq. 

TG\s = TG^^s- (3.21) 
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will be useful in the following discussion. 

The symmetrized bound state vertex function |r)s = ^ajF) = A^J^il^^) satisfies the 
equations 

|r), = -VG%s\n = -SAsV'G'ssin 

,(r| = -s{r\G%sV = -3 s{t\g'bsV'A3 , (3.22) 

where we used 

$:c^,,|r), = 3|r),. (3.23) 

i 

The normalization condition for this vertex function follows readily from the previous deriva- 
tion 

1 = .(r| -^ n - siriGls ^ Glsin = -s{rns ^ (vcy |r). . (3.24) 

Note that this equation has a structure identical to all of the previous results. 

To make a better comparison with the spectator equations, it is convenient to obtain 
the BS equations and normalization condition for the symmetrized subamplitudes. First 
introduce the operator 

A, = ^ (1 + CVjk) CViv = ^ CVw (1 + CVrk') . (3.25) 

The two forms of Aui are identical regardless of the specific values of i and i' . If z = i', then 
Vii' = C ^^^ 'Pjk = 'Pj'k' because both sets of indices {ijk} and {i'j'k'} must assume the 
values 1,2, and 3. li i y^ i', then i' must equal either j or k. For definiteness assume i' = j. 
Then we may use the identity 

VjkVi, = VijV.k (3.26) 

to prove the equivalence. An identical argument works if i' = k. These identities will be 
used frequently in the discussion which follows. Note that 

Y,Au' = A3 = Y.Au'. (3.27) 

i i' 

Now we define the symmetrized subamplitudes T** and |F*)s 

T^' = Y.AuT'''Ae, (3.28) 

\T% = J2Ae\r'). (3.29) 



Using ( |3.27| ) gives the relations 
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Y^T' = A.Y.'^'' A,, (3.30) 

Y. T"' = ^3 E ^'''^3 = A^TA^ = T (3.31) 

ii' w 

E r)s = ^3 E ir') = ^3|r) = |r), . (3.32) 

i I 



From the definitions it follows immediately that 



p.^T'' = (T'' = T^'Vj'k' (3.33) 

V,k\T% = C\T%. (3.34) 



It also follows from ( p.25|) that 



Vi,\T% = C\T^%. (3.35) 

To prove these relations define, for some operators O^, O"^, O^, the operator C 

^' = Y^^iO'. (3.36) 



Then 



v,,d' = ^v,,ii + cv,k)J2C^ii<^' 



^ (1 + (v^k) V,, (o' + cv^p' + CV^kO^ 

^ (1 + C,V^k) {ViP' + CO' + C,V^kV,kO^ 



6 



6 

CO^ . (3.37) 



In the last step we used {1 + CPik)'Pik = C (l + C^ife)- The first of the relations (|3.35|) 
follows immediately from (|3.37| ), and the second by taking the hermitian conjugate. 
Next, note that 



AuV' = Ay = V'Au, (3.38) 



and 






^V'E^*^^^*' = C^V%, (3.39) 



where in the last step we have used 



^(1 + CPJk)Y.^^iVu' = \{l + CPjk)KPu' ■ (3.40) 



Starting from Eqs. ( |3^ ) or (p.5|), using the definition ( p.28| ) and the properties ( |3.39| ), ( p. 381 ), 
( p.30| ), and (|3.27| ) gives the following equations for the symmetrized subamplitudes 

T-' = I C ^F'Grip,, - v'G'j^s E T'' (3.41) 

= I C \Vu'V'gt,^ - Y: r'G%sV' . (3.42) 



Summing Eq. (|3.41| ) over i and i' recovers Eq. (|3.18|) , demonstrating consistency. Note that 
Eq. ( |3.41D is equivalent to 

T" = i -V'Gt^ - V'G'bs E T^' (3-43) 

T"' = C T'Vii' , (3.44) 

where ( p.44| ) merely recovers the symmetry ( p.35[ ). Using this symmetry again, we may 
reduce ( p.43|) to an equation for the single amplitude T** 

T- = I ^-V'Gr^ - V'G'bs (1 + C^i, + C^ifc) T" 

= ^ \v'Gr^ - V'G'^s (1 + CV^, + CVjkV.,V,k) r' 

= ^ \v'Gr^ - V'G'ss (1 + 2C^.,) T' ■ (3.45) 

Hence it is sufficient to solve one equation for i = i', and obtain all other amplitudes from 
( p.44| ). Using Eq. (p.l8|) for the symmetrized two-body scattering amplitude, the equation 



for T*' may be further reduced to 

T" = i -M'Gr^ - 2C M'G^gVijT' . (3.46) 

Similarly, the equations for the symmetrized subvertex functions are 

r\ = -V'G^^s (1 + 2C^.,) \r% = -V'G^n (3.47) 

= ~2C M'G'j,s^,,\T%. (3.48) 

Using the definitions of the symmetric subvertex functions, the normalization condition 
( p.24| ) can be further simplified. First note that 



' yG%s)=A,4;,IVG%s)ECP.'--^A'^G-,, + v"^}j:CV.,. (3.49) 



To obtain this we used Eqs. (|3.2(]|J3.2T| ) and the fact that the exchange operators Vij are 
independent of P^, so that dVij/dP"^ = for any exchange operator. Substituting this 
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relation into (p.24| ) and using the bound state equation (|3.47|) and ^sir)^ = \T)s gives 
immediately 

1 = -3z ,(P| G, ^^ |r), + 3^ ,(r| G^^sG^ ^G'bs \n ■ (3-50) 

This condition can be expressed in terms of a single subvertex function, which we will 
choose to be ir*)^. Denoting the operator G^^gGi {d V /dP'^) G^^ by Og, and remembering 
that this operator contains a factor of 1 + (Vjk, the second term reduces to 

,(r| oi \r% = ,{r\ (i + (v,, + cv,k) oi (i + (v,, + (v.^) \r% 

= ,(P| (1 + 2CP.,) O^ (1 + 2CP.,) IF), (3.51) 

where we used VjkVijVjk = Vik -^ Vij as the operators Vjk are eliminated using ( p.34| ) and 
P|^ = 1. Similarly, since the operator Vjk commutes with 0\ = Gi{d G'^^g/dP'^), the first 
term becomes 

,(r| o\ \T% = ,{r\o\ (1 + cv,, + cv,k) \r% 

= s{r\o\{i + 2cn,)\r)s. (3.52) 

Combining these terms gives a normalization condition expressed in terms of |r*)s 

1 = -3ts{r\ G,^-^ (1 + 2QVi,) \T% 

+3z,{r\ (1 + 2Cn,) G'^sG.. -qPtGbs (1 + 2C^u) l^s ■ (3.53) 

We are now ready to turn our discussion to the three-body spectator equation. 



B. Normalization of the Three-body Gross Vertex Functions 

The three-body spectator or Gross equation for identical particles is obtained from the 
symmetrized three-body Bethe-Salpeter equation ( p.46| ) by placing the spectator particle on 
its positive energy mass shell ^^. In order to obtain a closed set of coupled equations a 
second particle must be placed on its mass shell in the initial and final states, as illustrated 
in Fig. p. Therefore, we label those particles that are on-shell by replacing subamplitudes 
T^i' by T**' ^ T]]: = Q^QjT'''Qi,Qj>. Noting that VjfQjVjf = Qy, and using Eqs. (^331) 
and ( |3.35| ) we have 



V^.T„ 



Vjkr;';, = QiQkV.kT'' Q, Qy = (Tl^, . (3.54) 



Similar relations exist for the initial state. Hence any amplitude can be obtained by 
the action of the permutation operators on the canonical amplitude Tj*, and it is necessary 
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+ 2C^ 




FIG. 9. Diagrammatic representation of the Gross equation for the amphtude Tj^g^. On-shell 
particles are labeled with an x . The final state particles emerge from the left side of each diagram, 
and putting the final state spectator on-shell consistently automatically forces a second particle in 
the final state to be on-shell. In this way two particles are always on-shell. 



only to find and solve the equation for this one amplitude. Note that we can choose i and 
j arbitrarily as long as i ^ j. The subamphtudes Tj*/ cannot be added unless the same 
pairs of particles are on shell in both the initial and final states. Therefore, there is no 
total scattering amplitude such as (p.7|) for the 3-body spectator equation. The underlying 
physical justification for the use of the spectator equations is described in detail in Refs. 
|0J^. We will show how to use the subamphtudes Tj*, in constructing the matrix elements 
in a forthcoming paper ||13|| . 

The equations for T**, were first derived in 0] in terms of the two-body scattering matrices 
M* •, by analysis of ladder and crossed ladder exchanges. Formally they follow from the three- 
body equation in the Bethe-Salpeter framework by appropriate restriction of the propagation 
of the on-shell particles. Thus, equations for Tjj can be obtained from Eq. ( ^.46 ) by replacing 
the two-body propagator G^^ by G*, where G* = GkQj- The inverse propagator in the 
inhomogeneous term, whose role was to cancel one of the spectator propagators when the 
equations are iterated, must be replaced by unity, iG^^ -^ 1, since the spectator is always 
on shell. This gives the following equations for Tj* 



T'i = -M]. - 2C Mi.G)Vi,T'i 



n 



n 



n 1 "^1 33 



Mi,-2CT'iVi,G]M: 



33 



33 ^3^j jj 



(3.55) 

(3.56) 



where the two body scattering operator M*- = QjM^Qj satisfies 



33 
33 



V 



33 



V G'M'- 

•^ 33 3 33 



V ■■ - M'G'V ■, 

33 33 3 33 ' 



(3.57) 
(3.58) 



Eq. ( |3.55| ) is illustrated diagrammatically in Fig. ^. In terms of the kernels V^ this equation 
becomes 



rpll 

33 



^33-T;}a + 2cv.3)G]V33- 



(3.59) 
(3.60) 
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Subvertex functions acquire an additional index to denote the second on-shell particle, 
|r*) — * |r*) = QiQj|r*), and the symmetrized subvertex functions satisfy the following 
equations 

\Vi) = -2CM]p)V,,\V)) (3.61) 

= -V]p]{l + 2CV,,)\T]). (3.62) 

Derivation of the normalization condition proceeds as for the Bethe-Salpeter case, but for 
clarity we will present the full derivation. Obtaining V ^^ from Eq. ( |3.6U| ), and substituting 
this into Eq. (|3.59| ) gives the following nonlinear equation for Tj* 



Tf, = ^],- ^f,G]S,irf^ - 3Tp,,G] V;^G)S,,T;] , (3.63) 

where Sij = (1 + 2C,'Pij). The existence of the bound state implies that 



in) (n, 



where, as before, i?*^- is regular at the pole at P^ = M^. Requiring the coefficient of the 
double pole on the rhs of this equation to be zero gives the bound state equation ( p. 62] ). 
Equating the coefficients of the single pole terms gives an equation similar to 



in) (HI = 3 lim in) 



(X]\{i+s.,G)v;^G)S.,\r)) 



(ni 



M2 - P2 

-I^})(^}|(l + 5.,G}F^)G55.,7^5^|, (3.65) 

where we used the fact that Sij commutes with QjG*. The bound state equation ensures 
that the terms involving TZJ, and any terms coming from the derivatives of the subvertex 
functions |rp with respect to P^ are zero. Since dSij/dP^ is also zero, the condition (p.65|) 
reduces to 

i = -3{(r}|5.,p|r5) + (r}|5.,^(G}T^q)5.,|r})| . (3.66) 

Distributing the derivative over the second term and using the bound state equation gives 
our final result for the normalization condition for the subamplitude |r*) 

1 = 3 |(r}| (1 + 2CP.,) (q)' iq) - (r}| (1 + 2C7^,) G] (y;^'G) (1 + 2C^,) |r})} . (3.67) 

As in the two body case, this equation could have been obtained directly from Eq. ( p.53| ). 

It is interesting to compare our principal result ( p.67| ) with the nonrelativistic normal- 
ization condition. If |\1/) = J2i 1"^*) is the total three body wave function, and we choose to 
express the normalization in terms of the i = 1 component, the normalization condition is 
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= 3(#| (1 + 2CP,,) m = 3(V^^| (1 + 2CP12) 1^') . 



(3.68) 



This is very similar to the normahzation condition for the Gross vertex function, 
Eq. ( |3.67| ). Adopting the convention that i = 1 and j = 2, using the fermion propaga- 
tor given in Eq. ( p.30|) , with k^ = P — ki — k2, gives (in the rest system) 



Glik,)] = G'.ik,) = G,{k; 



2P2 



G-sik; 



2M, 



^3(^:3) 7° ^3(^3) 



(3.69) 



B 



Neglecting the terms involving the derivative of the kernel in Eq. ( |3.67| ), decomposing the 
propagator into positive and negative energy parts, 



G'3(A;3) 



2Ei 



ks , 



E 



M(k3,A)M(k3,A) 



v{-k3,X)v{-k3,X) 



Eki + -Efca + ^ks — Mb -Efcs + Mb — E^^ — E, 



k2 



(3.70) 



keeping only the positive energy part and using u(k3, A')7°M(k3, A) = 2Ek^6x'x, allows us 
to reduce the normalization condition (|3.671 ) for the bound state vertex function for three 
identical fermions to 



Sj:{^l\{l + 2CVu)\^l) 



(3.71) 



where the wave function is related to the vertex function by 

1 u{k„X)\Tl) 



\^i) 



v/2M^ Ek, + Ek, + E, 



fes 



Mf 



(3.72) 



This correspondence will be developed in greater detail elsewhere. Here we only wish to 
emphasize that the covariant normalization condition ( |3.67|) reduces to the expected non- 
relativistic limit. 



IV. CONCLUSION 



The principal results of this paper are the derivation of the normalization conditions 
( p.24| ) and ( |3.53| ) for the three-body BS vertex functions, and ( p.67| ) for the three-body 
Gross vertex functions. These results are new. 

In the course of deriving these results we also showed that the normalization condition 
for the two-body Gross equation is identical to the requirement that the charge of a bound 
state be equal to sum of the charges of its constituents. Since we derived the normalization 
condition without any reference to the charge, this is equivalent to a proof that charge is 
automatically conserved by the spectator theory. This result has been used in previous work 
0, but the proof is new. A similar proof should be possible for the three body system, and 
will be presented in a subsequent paper [1T3| . 
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APPENDIX: FEYNMAN RULES 



In this paper, we use the Feynman rules as given in ||15|| . 

The variation from the usual field theoretical practice is with regard to the normalization 
of the symmetrized two- and three-body scattering matrices M and T. The practice in most 
field theory texts is to use an unnormalized symmetrization operator and to include all 
information about the symmetry of the scattering process in an overall phase space factor. 
Here we choose to use normalized symmetrization operators A2 for the two-body case and 



^3 for the three-body case, as defined by (|3.16|) . This results in integral equations for the 



scattering amplitudes which are consistently of the form of the nonrelativistic Lippmann- 
Schwinger equations and preserves the nonrelativistic normalization of the wave functions. 
Any matrix element that involves an occurrence of the symmetrized scattering matrices as 
used here can be used in the usual expression for cross sections by making the substitutions 
M ^ 2M and T ^ 6T. 
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FIGURE CAPTIONS 



Figure 1: Diagrammatic representation of the two-body BS equation for the scattering 
matrix. 



Figure 2: Diagrammatic representation of the symmetrized kernel. 

Figure 3: Diagrammatic representation of the two-body Gross equation for the scattering 
matrix. The x on the hne for particle 1 indicates that it is on shell. 

Figure 4: Diagrammatic representation of the symmetrized kernel for Gross equation. 
As in Fig. 3, the x indicates that the particle is on shell. 

Figure 5: Diagrammatic representation of the normalization condition for the two-body 
Gross vertex function. The dashed line represents the derivative d/dP'^ , and the x means 
that the particle is on shell. 



Figure 6: Diagrams which contribute to the electromagnetic charge of the two body 
bound state. 



Figure 7: Location in the pio complex plane of the 6 poles from the three nucleon 
propagators in Eq. (2.34). The last two terms in Eq. (2.33) emerge if the integral of pio over 
the contour C enclosing the poles 3a and la is evaluated using the residue theorem. 

Figure 8: Diagrammatic representation of the Faddeev equations for the amplitudes T^*. 
Note that the spectator is identified by the solid dot. 

Figure 9: Diagrammatic representation of the Gross equation for the amplitude T22- 
On-shell particles are labeled with an x . The final state particles emerge from the left side 
of each diagram, and putting the final state spectator on-shell consistently automatically 
forces a second particle in the final state to be on-shell. In this way two particles are always 
on-shell. 
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